Power law load dependence of atomic friction 



o 
o 

u 

Oh- 
< 
(N 



O 



-a 
c 

o 
o 



> 

o 
o 
o 

-I— > 



i 

C 

O 

o 



C. Fusco* and A. Fasolino 
Theoretical Physics, NSRIM, University of Nijmegen, 
Toernooiveld 1, 6525 ED Nijmegen, The Netherlands 
(February 2, 2008) 

We present a theoretical study of the dynamics of a tip scanning a graphite surface as a func- 
tion of the applied load. From the analysis of the lateral forces, we extract the friction force and 
the corrugation of the effective tip-surface interaction potential. We find both the friction force 
and potential amplitude to have a power law dependence on applied load with exponent ~ 1.6. 
We interpret these results as characteristic of sharp undeformable tips in contrast to the case of 
macroscopic and elastic microscopic contacts. 



It is well known that macroscopic friction is propor- 
tional to the applied load but the load dependence of 
atomistic friction is still under investigation and has been 
discussed in a limited number of experimental and the- 
oretical works [1-8]. Usually the load dependence is de- 
scribed by means of contact mechanics continuum theo- 
ries [7,8] which do not give information about the atomic 
interactions in the sliding contact region and assume a 
spherical tip. Moreover, it has been argued that, de- 
pending on the shape of the tip, power laws with differ- 
ent exponents can be found [7] . In this Letter, we present 
a detailed study of the load dependence of atomic-scale 
friction in the case of a sharp tip-surface contact, find- 
ing a power law dependence with exponent larger than 
one. We also show how the effective tip-surface potential 
energy barriers can be derived from force measurements, 
thus establishing a connection between friction and the 
interatomic potential. Recent works using noncontact 
mode Atomic Force Microscopy (AFM) have shown the 
possibility to reconstruct the tip-surface potential [9] , but 
they lack the link between the corrugation of the poten- 
tial and the friction force. The ideal way to achieve this 
goal is provided by a study of the load dependence. 

We address the load dependence of friction by simu- 
lating the dynamics of a tip scanning a rigid monolayer 
graphite surface, extending to three dimensions (3-D) the 
Tomlinson-like models of AFM [10]. Experimental evi- 
dences show that the tip usually cleaves small graphite 
flakes which remain attached to it [11]. Since the con- 
tact diameter between the tip and the surface can be 
very small we consider the limiting case in which the tip 
is formed by a single carbon atom connected via har- 
monic interactions with force constants K x , K y and K z 
to a support moving along the scanning direction. The 
carbon atom of coordinates (x, y, z) interacts with the 
graphite surface via the specific empirical potential Vrs 
for graphite given by [12]: 



where 9(r) is the Heaviside function, rxj is the distance 
between the tip carbon atom and the j— th substrate 
atom and 



Vi(r) = e t (V2ft(r-r ) _ 2e -0 i{r -r o )\ + y . 



1,2 



(2) 



with vi = ei — 62 and v 2 — 0. The values of the param- 
eters are r = 0.371 nm, ei = 5.355 meV, e 2 = 2.614 
meV, fa = 14.693 nm" 1 and (3 2 = 21.029 nm -1 . The tip 
support is moved with constant velocity v scan along the 
scanning line x scan = v scan t, y scan =constant. We solve 
numerically the equations of motion in the constant force 
mode, i.e. we set K z — and add a constant force Fi oa( i 
in the downward z direction, including also a damping 
term proportional to the atom velocity, which takes into 
account dissipation mechanisms: 
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We adopt an atomistic approach, assuming for m the 
mass of a single carbon atom (m = 1.92 • 10 -26 kg) and 
for the damping parameter 77 = 1 ps _1 , which is a value 
appropriate for dissipation of energy and momentum at 
the atomic scale (see e.g. [13]). Here we show results for 
K x — K y = 1 N/m which are typical values of AFM, 
whereas our scanning velocity v scan = 0.4 m/s is much 
higher than in experiments. Our choice of parameters 
makes the dynamics underdamped. 
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FIG. 1. Tip-surface potential Vrs as a function of the 
tip-surface distance z with x and y coordinates taken as those 
of the hollow site (solid line). Points give the values of Vrs at 
several positions along the actual trajectory as a function of 
the corresponding instantaneous value of z. Results for all in- 
dicated loads are obtained for simulations where the scanning 
direction is along a row of atoms (y S can = 0.17 nm). 
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FIG. 2. Contour plot of the potential Vrs for Fi oa d = 1 
nN (top) and Fi oa d = 4 nN (bottom) with z fixed at the av- 
erage value < z > of the height obtained by the simulations 
(< z >= 0.25 nm for Fi oad = 1 nN and < z >= 0.2 nm for 
Fioad = 4 nN). The minimum V m in and the maximum V m ax 
of Vrs are located respectively at the hollow site and on top 
of one atom. We show 10 contour lines between V m in and 
Vmax, separated by an energy interval A. V m in = 123.5 meV, 
Vmax = 158.5 meV and A = 3.2 meV for F toa d = 1 nN, while 
Vmin = 774 meV, V max = 1.01 eV and A = 23.7 meV for 
Fioad = 4 nN. The horizontal thick solid line indicates the 
scanning direction (y 3ca n = 0.17 nm), while the crosses are 
points along the actual atomic trajectory taken every 4 ps. 



The energetics underlying the tip motion for different 
loads is illustrated in Figs. 1 and 2. In Fig. 1 we present 
the tip-surface potential Vrs as a function of the verti- 
cal tip-surface distance z, for in-plane coordinates cor- 
responding to the position of the hollow site, on which 
we superimpose the values of Vrs at actual positions of 
the tip during the simulations for different values of the 
load. The tip position probes the attractive part of the 
potential for low loads (Fi oa d = 0.1 nN) and moves closer 
and closer to the repulsive core of the substrate atoms for 
increasing loads. Clearly the whole tip-surface potential 
is probed by the motion of the tip for different loads, as 
suggested in Ref. [2]. In Fig. 2 we show the contour plots 
of Vrs , calculated for z given by the average value of 
the tip height at two values of the load. The scanning 
line Dscan = 0.17 nm corresponds to scanning along a 
row of atoms. The actual trajectory is also reported: for 
the smaller load the motion follows a zig-zag-like pattern, 
while the trajectory acquires a one-dimensional stick-slip 
character for the higher load (see also [2]). 



Fig. 3 presents typical force plots as a function of x scan 
and x [14] for Fi oad = 2 nN and Fi oad = 4 nN. The saw- 
tooth behavior characteristic of stick-slip motion is deter- 
mined by the competition between the elastic force F e ; 
and the force Fts due to the tip-surface potential Vrs- 
Elastic energy, accumulated in the spring, is counterbal- 
anced by the substrate attraction, until, suddenly, the 
tip slips to another minimum. Therefore, while sticking, 
F e i mirrors Fts- This fact can be used to derive Vrs 
itself. The solid lines in Fig. 3 represent the lateral force 
along the scanning direction, as obtained by our sim- 
ulations. Increasing the load enhances the sawtooth be- 
havior of the stick-slip motion and gives rise to a larger 
initial sticking region, often observed experimentally (see 
e.g. [15]). 

As shown in Fig. 2 the actual trajectory does not neces- 
sarily follow the scanning line so that the potential energy 
landscape during the motion is not known a priori. How- 
ever, we can extract the effective value of the energetic 
barrier Vo for a given Fi oad directly from the force plots. 
By assuming a sinusoidal shape of Fts [5,16] and noting 
that Fts should average to zero for a periodic substrate, 
we can reconstruct Ff s . 
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FIG. 3. Lateral forces along the scanning direction for 
Fi oad = 2 nN [(a),(c)] and Fi oad = 4 nN [(b), (d)], plotted as a 
function of x scan [(a),(b)] and x [(c), (d)], for y SC an = 0.17 nm. 
Solid lines are the elastic forces F^ obtained by simulations, 
while the dashed lines in (c) and (d) are static calculations of 
the tip-surface force (with reverted sign) —F^ s at (y, z) de- 
termined by averaging y(t) and z(t) given by the dynamics. 
The dotted lines in (a), (b) give the slope K e ff of the sticking 
part (K eff = 0.78 N/m for F load = 2 nN and K eff = 0.89 
N/m for F load = 4 nN). 



Then, the tip-surface potential Vts is simply given, up 
to a constant, by integrating F^ s : 



V TS {x) 



j F^ s dx 



V . /2ttx\ 

— sin 

2 \ a J 



const. (4) 



with Vo = Foa/{2n), where a and Fo are the period and 
the amplitude of Fts- In Figs. 3(c)-(d) we show —F^ s 
(dashed lines) obtained by static calculations for (y, z) 
determined by averaging y(t) and z(t) given by the dy- 
namics. Indeed, — F^, s follows the sticking parts of F^ 
quite well. The first stick signal, of larger amplitude, is 
the most suitable to estimate the amplitude of F^ s . 
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FIG. 4. Load dependence of the energy barrier (a) and of 
the friction force (b) for y sca n = 0.17 nm. Inset of (b): re- 
lation between friction force and energy barrier. The solid 
circles are the result of the simulations and the dashed lines 
are power law fits, as explained in the text. 

The resulting dependence of the energy barrier Vb on 
the load is shown in Fig. 4(a). An increase of Vo with 
Fioad has also been found experimentally [2,17]. The 
dashed line in Fig. 4(a) is a power law fit to the numer- 
ical data with exponent <~ 1.6. Fig. 4(b) illustrates the 
lateral friction force Ff r i C , obtained as the mean value 
of the instantaneous lateral force F^, as a function of 
Fioad- Also these data can be fitted by a power law with 
exponent ~ 1.6. Thus, the linear relation between macro- 
scopic friction and load does not hold at the microscopic 
level. Moreover, the exponent is different from the 2/3 
expected for a Hertzian contact [7,8]. As pointed out in 
Ref. [7] even a small deviation from the spherical shape 
of the tip can be responsible for a change in the power 
law exponent. We conjecture that exponents larger than 
one are characteristic of sharp, undcformable tip-surface 
contacts. Note that, since the power law exponents for Vo 
and Ff r i c as a function of Fi oa d are very close, an approx- 
imately linear relation between F/ r j c and Vo is recovered 
(see inset of Fig. 4(b)), indicating a direct link between 
atomistic friction and energy barriers for diffusion. 

In conclusion, we have presented a theoretical study 
of the 3D dynamics of a tip scanning a graphite surface 
with realistic interactions as a function of the applied 
load. We predict a power law behavior with exponent 
~ 1.6 for the friction force as a function of applied load, 
at variance with macroscopic behavior. The study of the 
load dependence establishes a direct linear relation be- 
tween friction and potential corrugation in the contact 
region. 
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